Stochastic Schrodinger equations with coloured noise 
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A natural non-Markovian extension of the theory of white noise quantum trajectories is pre- 
sented. In order to introduce memory effects in the formalism an Ornstein-Uhlenbeck coloured 
noise is considered as the output driving process. Under certain conditions a random Hamiltonian 
evolution is recovered. Moreover, non-Markovian stochastic Schrodinger equations which unravel 
non-Markovian master equations are derived. 
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Recently, stochastic wave function methods for the de- 
scription of open quantum systems have received consid- 
erable attention [2, H, HI]. These approaches have been 
mainly motivated by the continuous measurement de- 
scription of detection schemes in quantum optics, namely 
direct photo detection, homodyne and heterodyne detec- 
tion. In general, selective, indirect and continuous quan- 
tum measurements allow for a description of the open 
quantum system in terms of the stochastic evolution of 
its own wave function. Typically, the random trajecto- 
ries of the state vector involve quantum jumps or diffu- 
sion processes and are described by stochastic differential 
equations for the wave function \ip t ). 

The relationship to the more traditional approach to 
open quantum systems in terms of Master equations for 
the density matrix p t is easily established by realizing 
that the latter can be expressed as p t = K[\tjjt)(tpt\], where 
E denotes the ensemble average over realizations of the 
stochastic process \ipt)- Thus, the Master equation evo- 
lution can be reproduced by generating a large number 
of trajectories of the state vector. This unravelling pro- 
cedure, called Monte Carlo wave function method, has 
gained considerable importance for the numerical simu- 
lation of complex open systems Q- 

The procedure described above has been a major 
breakthrough in the description of the Markovian dynam- 
ics of open quantum systems, where it has found many 
applications. At present an active line of research is to 
find generalizations of the stochastic Schrodinger equa- 
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tion in the non-Markovian regime, which have a sound 
physical interpretation. 

Essentially, two strategies have been investigated to 
accomplish this goal. A first strategy is to start from 
non-Markovian master equations and to try to con- 
struct pure state unravellings. It is well-know that non- 
Markovian master equations can be constructed with the 
help of projection operator techniques. For example, 
time-convolutionless (TCL) projection techniques yield 
master equations that are local in time and allow for suit- 
able unravellings (0, 0] for details). In the same spirit, 
non-Markovian generalizations of the Lindblad theory Q 
have been proposed and shown to have corresponding 
non-Markovian unravellings 0, 0] • The measurement in- 
terpretation of these unravellings is still an open prob- 
lem. However, special cases of unravellings with a clear 
measurement interpretation have been shown to exist for 
discrete time models of repeated interaction [8[. 

The second strategy is first to generalize directly the 
stochastic Schrodinger equation and secondly to show 
if it provides an unravelling of a suitable Master equa- 
tion and if it has a measurement interpretation. In 
0) B B EH EH E3] starting from a measurement scheme 
a rigourous non-Markovian random evolution has been 
derived. The question of the possibility to express this 
random evolution in terms of pure states in the spirit 
of unravelling is still highly debated [t| Ell- For the 
usual scheme of indirect quantum measurement it has 
been shown that in general such an interpretation is not 
accessible [§, [lfl. Other investigations promise to get 
unravelling EH • 

In classical statistical physics the most straightforward 
way to introduce memory effects is to consider coloured 
noise [14] . This approach has already been successfully 
transferred to the theory of isolated quantum systems, 
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by introducing; coloured noise in the Hamiltonian of the 
evolution [3, H!, EH, G3]- I n this paper we show how to 
introduce memory effects in the stochastic Schrodinger 
equation with the help of the addition of coloured noise. 
Specifically, we will illustrate the approach by using the 
Ornstcin-Uhlcnbeck process. 

The alternative strategy, proposed here, is focused on 
an extension of the canonical Markovian approach. To be 
precise, generic Markovian stochastic Schrodinger equa- 
tions are justified by the use of a linear equation for a 
non normalized state ip t of the following form 



dipt = Aij t dt + BiptdX u 



(1) 



where the formal derivative of (X t ) is a white noise, that 
is, E[X t X a ] = S(t - s) and E[X t ] = (E denotes the 
expectation under a reference probability). In general, 
the process (X t ) can be a Brownian motion, a counting 
process or a mixing of both noises. Physically, the noise 
terms represent the output signal of the measurement 
and the probability of an event is traditionally described 
by ||V't|| 2 - Usually, asking \\tpt\\ 2 to represent a prob- 
ability density (giving the physical probability) implies 
-^[llV'tll 2 ] = 1 an( i this last equality imposes the form of 
the operators A and B. Finally, a new stochastic differen- 
tial equation for normalized vector states {ipt) is derived. 
In the equation for the normalized states a new white 
noise appears, obtained by a suitable Girsanov trans- 
formation of the reference Brownian motion. A natu- 
ral introduction of non-Markovian memory effects in the 
framework sketched above, can be obtained by replacing 
the white noise (X t ) by a coloured Ornstein-Uhlenbeck 
noise. 

To keep things simple here we consider a one- 
dimensional Ornstein-Uhlenbeck process defined by X t — 
f o * e -T(*- s )dW s , where (W t ) is a one dimensional Brow- 
nian motion, under a reference probability Q. The 
Ornstein-Uhlenbeck process (X t ) satisfies the stochastic 
differential equation 



dX t = -~fX t dt + dW t 



(2) 



where 7 > 0. The high non-Mar kovianity of the process 
is emphasized by the fact that the correlation function is 
not a pure delta function, but we have EQ[X t X s ] = 5{t — 
s ) — 2 (e _7 '* _s ' + e _7 (* +s ') and the Markovian regime is 
recovered in the limit 7 [ 0. 

In the following, first we will show how this approach 
allows to reconstruct dynamical models with random 
Hamiltonians. Next, we will demonstrate how to con- 
struct a class of non-Markovian stochastic Schrodinger 
equations, that has an interpretation in terms of quan- 
tum measurement and unravelling. 

Random Hamiltonians. Actually, the random Hamil- 
tonian models arise naturally by translating the action 
of the Ornstein-Uhlenbeck process © in the definition 
of the linear stochastic Schrodinger equation JT]). With 
the output driving process (X t ) satisfying (QJ), Eq. (fT|) 



can be rewritten in the following form 

dyj t = (A - 1 X t B)ij H dt + B^ t dW t . 



(3) 



where the initial condition is a wave function ipo, such 
that \\i, \\ 2 = 1. 

According to the laws of quantum measurement the 
probability of the output is described by the density 
H^tll 2 . This yields a measurement interpretation of Eq. 
([3|). The quantity HV'tll 2 defines a probability density, if 
EqUI^D 2 ] = 1, is satisfied for all times t. The consistency 
of these densities needs also the process (HV'tll 2 ) to be a 
martingale (we refer to [X9\ for complete details). By ltd 
calculus rules, the stochastic differential of ||V't|| 2 turns 
out to be 



d(^|V*) = (dip t \i>t) 
= (i> t \[A*+A- 



+ (iPt\diPt) + (d^ t |d^ t ) 
■jX t (B* + B) + B*B] x/, t )dt 

+ {4> t \(B* + B)i> t )dW t . (4) 



In order to ensure that Eq^^I) 2 ] = I, for all times t, 
we must impose that the term in dt must be equal to 
zero (and this guarantees also the martingale property). 
Indeed, in this way, we get Eq^^II 2 ] = EqHIV'oII 2 ] = E 
since we have E^fdWi] = 0. It follows that we must have 



.4* 



A - jX t (B* +B) + B*B = 0, 



Vi. 



(5) 



This imposes that there are two self-adjoint operators K 
and H such that B = — iK and A = — iH 
Eq. becomes 



±K 2 . Then, 



dipt 



1 



-i (H - ~/X t K) 



iPtdt-iKiPtdWt. (6) 



Let T^ denotes the time ordering exponential, the for- 
mal solution of Eq. ^ is given by 



i / (H - 1 X S K) ds-i KdW s \ V, 



In other terms, the evolution of the quantum system is 
completely determined by the time-dependent, random 

Hamiltonian H t = H + (\V t - jX t ^j K (a formal expres- 
sion, due to the presence of W t ). Note that Eq. ^ gives 
d|iV't|| 2 = 0) m agreement with a purely Hamiltonian 
evolution. This shows, that the usual measurement in- 
terpretation of |T]) coloured with an Ornstein-Uhlenbeck 
process gives raise to a random Hamiltonian evolution. 
The property ||V't|| 2 = 1 implies that we have not ex- 
tracted information according to our measurement inter- 
pretation. We recover the framework of the evolution 
of an isolated closed system incorporating a random en- 
vironment characterized in terms of Ornstein-Uhlenbeck 
noise. 

In order to investigate the property of Eq. ([6]) in terms 
of unravelling, it is interesting to consider the evolution of 
the corresponding density matrices. To this end, we con- 
sider the pure state process (pt) defined by pt = \tpt){tpt\- 
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We define also the mean state rj t = Eqlpt]. Using Ito 
rules, the process (p t ) satisfies the stochastic differential 
equation (SDE) 



dpt 



-i[H - jX t K, Pt ]dt - i[K, pt}dW t 



~[K,[K,p t ]]dt, (7) 



which is, of course, equivalent to ©. Let us stress that 
the presence of the Ornstein-Uhlenbeck process implies 
that the solution (p t ) of Eq. ((7]) is not a Markov process. 
For the mean state (j] t ) we have then 

^r, t = -i[H, m ] - \ [K, + i 7 [if, E Q LY t(0t ]]. (8) 

The above equation can be naturally considered as a Mas- 
ter equation, but its particularity is that it is not a closed 
equation for the mean state r\t ■ Actually, we have derived 
a model with memory for the mean state. Indeed, the 
term i7 \K, Kq[Xtpt]\ introduces non-Markovian memory 
effects in the dynamics. Moreover, Eq. ([6]) is an unravel- 
ling of the Master equation ([5]). 

Random coefficient The random Hamiltonian case we 
discussed above arose only because the coefficients A and 
B in the Eq. |T]) were assumed to be not random. In the 
following we will address the situation, in which A and B 
are assumed to be functions of the Ornstein-Uhlenbeck 
process (X t ), thus becoming random operators them- 
selves. In other words, we consider now a SDE of the 
more general form 



d^ t = A{X t )i> t dt + B(X t )i> t dX t . 



(9) 



Inserting the definition ([2]) of the Ornstein-Uhlenbeck 
process into Eq. §§§ yields 

d^ t = (A(X t ) - 7 X t B{X t ))^ t dt + B(X t )i> t dW t . (10) 

Again, asking the norm ||V't|| 2 to be a probability density 
the condition ([5]) is translated into the following one 

A*(X t ) + A(X t ) - 1 X t (B*{X t ) + B{X t )) 

+ B*(X t )B(X t )=0. (11) 

In contrast to the previous case with deterministic co- 
efficients A and B, we can now establish a relationship 
between the operators A(X t ) and B(X t ). More precisely, 
the condition (fTTj) implies that there exists a self-adjoint 
operator H(X t ) such that 

A(X t ) - jX t B{X t ) = -iH(X t ) - i B*(X t )B(X t ). 

Then, Eq. HU]) becomes 



d^ t = ( -iH(X t ) - - B*(X t )B(X t ) ) ^ t dt 



In the general case d||V't|| ^ and the above evolution 
cannot be interpreted in terms of a random Hamiltonian. 
This is the signature that the above evolution describes 
the extraction of information from the system in terms 
of indirect quantum measurement. Indeed, the measure- 
ment interpretation relies upon the property that the 
norm ||V'|| 2 defines a probability density. According to 
this property the relevant physical probability is defined 
by 

P T (d^) = U T (uj)\\ 2 ^{dLa), 

where Q is the underlying reference probability as intro- 
duced in the definition of the Ornstein-Uhlenbeck pro- 
cess. In other words, the probability P T is the physical 
probability law of the events which could occur in the 
time interval [0,T]. When the time T varies, the proba- 
bilities P T are "consistent", because HV'tII 2 is a "martin- 
gale" [HI . Consistency means that if we take < S < T 
and an event determined by conditions only in the time 
interval (0,5), then P T (F) = ¥ S {F). Another impor- 
tant property of these probabilities is that they can be 
expressed in terms of quantum expectations on the pre- 
measurement state of positive operator valued measures, 
as stated by the axioms of quantum mechanics (l9j . 

In order to describe the stochastic dynamics of the sys- 
tem undergoing continuous indirect measurement we de- 
rive the corresponding stochastic Schrodinger equation 
for the normalized state vector ifj t = V't/IIV'tll- Usually, 
such a process is called a quantum trajectory and satis- 
fies a non-linear stochastic differential equation. To de- 
rive this equation we need to describe the driving random 
process in terms of the new physical probability P T . The 
famous Girsanov theorem gives that, under P T , the pro- 
cess 



W t = W t 



m(s)ds, 



with 



m(t) = ($ t \(B*(X t )+B(X t ))$t), 



is a new Brownian motion. This allows to express the 
output of the measurement X t (which was an Ornstein- 
Uhlenbeck coloured noise under the reference probability 
Q) in the form 

X t = [ e'^ t " s '>m{s)ds+ [ e^*-^dW a . 
Jo Jo 

Moreover, by using Ito calculus we can write the non- 
linear stochastic differential equation for the normalized 
state vector -0 as 



d<<p t 



B{X t ) - 
m(t) 



m(t) 



iptdW t 



iH(X t ) 



B{X t )^tdW t . (12) 



B(X t ) + -B*(X t )B(X t ) 



m(t) 2 



iptdt. (13) 
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Let us stress that the presence of the coloured noise X t in 
the above equation makes the process ip non-Markovian. 
Thus, Eq. (fT3)) is an example of a non-Markovian stochas- 
tic Schrodinger equation that admits a clear interpreta- 
tion in terms of indirect quantum measurement. Eq. (fl~3|) 
has the usual structure of a stochastic Schrodinger equa- 
tion. Its particularity is the fact that the operators H and 
B are now random functions of the Ornstein-Uhlenbeck 
process X t . 

We shall investigate the evolution of the corresponding 
density matrix. To this end we define 

Pt = \i>t)(i>t\- 

The non-linear stochastic master equation for the density 
matrix p t takes the form 

Apt = C(X t )[pt]dt + [B(X t )pt + p t B*(X t ) 

- Tr{[B(X t ) + B* [X t )] p t } Pt] dW t , (14) 



where 



= -i[H(X t ),p] 

-±{B*(Xt)B(Xt),p} + B(Xt)pB*(Xt), 



defines a random Lindblad generator. Note, that again 
the originality relies on the fact that the operators H and 
B are random (see [H, for similar considerations). 

The mean state r\ t = ^r T [pt] = ^q[\ijjt) {tpt]] satisfies 
the following non-closed equation 



-77 t =E P T [C(X t )[pt 



(15) 



This equation describes a non-Markovian deterministic 
Master equation. Of course, by construction the pro- 
cess -ipt yields a non-Markovian unravelling of the Master 
equation (fT5|) . 

A special case is featured by considering B to be con- 
stant. In this case, we recover a model for the open sys- 
tem with a random Hamiltonian and a usual dissipative 
part in Lindblad form with a non random operator B. 

In conclusion, a special class of non-Markovian 
stochastic Schrodinger equations and stochastic master 
equations has been derived. The framework of the tradi- 
tional Markovian white noise theory has been extended 
by introducing the use of coloured Ornstein Uhlenbeck 
noise. This natural extension allows to describe random 
dynamics for open quantum system with memory effects. 
Two important categories of non-Markovian models has 
been derived. On the one hand, in a natural way, we have 
recovered the special case of random Hamiltonians. On 
the other hand, a generalization of Lindblad type equa- 
tions has been obtained. The memory effects, influenced 
by the Ornstein Uhlenbeck noise, are encoded in the def- 
inition of the random operators which define either the 
Hamiltonian or the Lindblad operator. Our approach, 
based on the description of stochastic Schrodinger equa- 
tions in terms of pure states, give rises to non-Markovian 
unravellings of these models. 
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